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Abstract. The hedgehog Skyrme model on three-sphere admits very rich spectrum of solitonic solutions 
which can be encompassed by a strikingly simple scheme. The main result of this paper is the statement 
of the tripartite structure of solutions of the model and the discovery in what configurations these solutions 
appear. The model contains features of more complicated models in General Relativity and as such can give 
insight into them. 
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1. Introduction to the Skyrme model on S 3 and motivation 

In this paper I discuss static solutions of the Skyrme model on the three-sphere of radius I? as a physical 
space. In a sense, to be explained later this model includes the Skyrme model on flat space R 3 in the limit 
R —> oo. The model contains two scales of length which constitute a free dimensionless parameter whose 
value is crucial to the existence of solutions of a given type. There are known several models in which similar 
behaviour is observed ( e.g. jsj) and because of this it deserves investigation. The Skyrme model on S 3 is the 
simplest nonlinear realization of such models which admits very rich structure of solitonic solutions. 

The Skyrme model on flat space was originally introduced as minimal nonlinear realization of an effective 
meson field theory where low energy baryons emerged as solitonic fields. Then nucleons were introduced as 
topological excitations of these fields. The basic chiral field in that model is SU (2) valued scalar field U ( x ) 
which is related to the singlet meson field a and the triplet pion field 7 ? by 

U(x) = ^-(<r(x)+ia ( a ) 2 + (if ) 2 = f 2 . 

' TT 

The model is defined by the Lagrange density 

(1-1) £ = \fl Tr (L,L») + J-^Tr ([L^, L v ] [L M , L v ]), 

where 

Lfj, = U+d p U 

is su( 2) valued topological current. The model is characterized by two constants f n - pion decay constant 
and e determining the strength of the fourth order term which was introduced by Skyrme in [Q in order to 
ensure existence of solitonic solutions, but it can be introduced in quite different way - on geometrical level 
which I shall sketch now. 

The Lagrange density 0) has a nice geometrical origin. An m-parameter Lie group Q valued scalar field 
U is a map M 3 p —> U(p) £ Q. If x = x{p) is a local coordinate system on M and £ are co-ordinates of the 
Lie group about the identity of the group, then the field induces a map R™ 9 x —> £( 2 ;) £ R m and respective 
mapping of the tangent bundle TA4 to the Lie algebra T e Q - the tangent space to the Lie group Q at the 
identity. Each one-parameter subgroup U(t) of Q such that U{t\ + £ 2 ) = U{t\)U(t 2 ) and 17(0) = e defines 
a smooth curve U(t) with tangent vector U(t) = U(t)u £ T V ^Q where u £ T e Q is tangent to e, so that 
C/ —1 C7 (t) is an element of T e Q. This fact may be used to define the repere mobile of the Lie algebra at the 
identity by 

a*(£) := U-'diUig). 

By that, with each v £ T p A4 we can associate an element of the Lie algebra 

T P M 3v^ ai ({( x ))^- v »(x) £ T e g, 
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which defines a Lie algebra valued one-form Q = —U~ 1 d fJ ,U daT on A4 which simultaneously is a matrix of 
connection one-forms whose curvature identically vanishes dfl — A Q = 0. From the connection we construct 
secondary quantities 0 := dfl and, since A4 is a Riemannian manifold, *12 and *0 (where * is the Hodge 
operator). Now we can build invariants being full scalar products (*f2, fl) and (*0,0), taking the first as 
the sigma term and the second as the self-interaction term for our field such that, imposing on /i, A, and U 
physical dimensions, we arrive at 0 starting from 

J C(U , dU)V^g d 4 z := J pTr (*ft A ft) + ATr (*0 A 0 ). 

M M 

This equation defines, in generally covariant manner, the lagrange density C and generalizes the flat Skyrme 
model to an arbitrary curved space-time. Another point of view, though only for static field configurations, 
was presented by Krush in || on the level of density of energy, and can be recast into the following scheme. 
Taking, as a physical quantity, the strain tensor D ik = — ^Tr (U~ 1 diUU~ 1 dkU) being the Killing form for 
the Lie algebra of the compact matrix group SU( 2), we use invariants originating from the secular equation 
det(l + t,D) = Ckt k , which serves as generating polynomial in t. The invariants are respectively c\ = Tr D, 
C 2 = ^(Tr D ) 2 — ^Tr D 2 , which in particular reproduces the sigma and self-interaction terms in ([b|), and 
additionally C 3 = det(D) which is proportional to the density of topological charge. This approach generalizes 
the energy density for static solutions to an arbitrary curved three-geometry. (Aside remark: such invariants 
can be easily derived if one notices that det(A) = exp(Tr log(A)) what is obvious for diagonalizable matrices 
and, by virtue of continuity, must hold in general. Putting A = E + tD , expanding in the variable t, and 
remembering that for matrices D of rank n the Taylor series expansion in t must terminate at t n , we get all 
invariants for D). 

The lagrange density 0 ) leads to the action for a field U: 


S[U}= f C(U,dU)^jd 4 x. 
Jm 


! M 

Here we immerse the Skyrme model within the simplest possible curved space-time which, from now on, 
we assume to have the structure of the Cartesian product M. = R x S\ (S^ is the three-sphere of radius R) 
which is mapped into a subset of R 4 by the map R x 5^ 9 p —> T(p) = In this map the line 

element takes the form 

^2 _ ij .2 d 2 j „/,2 d 2 „:„2 I iq 2 , „:„ 2 „ qj ,„ 2 \ 


ds 2 = dt 2 - R 2 dip 2 - R 2 sin 2 ip {Ad 2 + sin 2 ddip 2 


and, disregarding singularities of the map, we assume that t £ (—00, +00), ip £ [0,7r], 1 ? £ [0, 7 r] and ip £ [0, 2n). 
For that space-time is invariant with respect to the group of translations in time, the first theorem of Noether 
assures for solutions the existence of a constant of motion - the energy integral. Moreover solutions which 
are static, and we limit ourself in this paper considering only such one, are critical points of the integral 


(9S3 ~w 

( 1 . 2 ) 


9 ik ) 


E = R d 


Js 3 


Too^/gs 3 ^ 3 ^ 


where Too is 00 component of the metrical energy-momentum tensor, which for static solutions takes the form 
(then Too = -£) 


(1.3) 


Too = 


1 ft 


1 


T r + lfie 2 /?' 1 (J s' i( 9s 3 ^ r (TiLj[L k , Li\). 


The integral 1.2 measures the energy of a map from the configuration space S 3 to the SU( 2) group which, as 


a Riemannian manifold, is equivalent to S 3 . So that we examine energies for mappings from S 3 into S 3 . 

The Euler-Lagrange equations for the Skyrme model on S 3 , defined above by the energy integral (0), 
lead to a system of three nonlinear partial differential equations of second order for thre e un known functions. 
To bypass this outstanding example of human endeavour, exploiting the invariance of (|l.2| ) under the action 
of SO(3) group of rotations, we restrict ourselves by considering only those fields which are S’0(3) invariant. 
According to the theorem given by Faddeev in || an extremum of a functional in the domain of invariant 
fields under the group of invariance of the functional is the extremum of the functional in the domain of all 
fields. The hedgehog Ansatz provides us with a subgroup of fields which are SO( 3) invariant. The hedgehog 
field, which in original Skyrme model describes baryons, is defined by 

(1.4) U{ip, ip, -d) = e i*°iHv,*)FW' 

and in language of Skyrme the pion field 7? is identified with the direction n(ip, i9) in isospin space and F(ip) 
measures its amplitude. The form of this Ansatz, of course, comes from the canonical parameterization of 
SU (2) group with addition of implicitly assumed equivariance of the map S 3 —> S 3 i.e. (ip, <p) —> ('I', 0, <!>) 
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such that >]) = F(ip), 0 = $, $ = (p. This is in accordance with assumptions made in 0 and by that we get a 
natural field-theoretical generalization of the model of harmonic maps between three-spheres of Bizon. Our 
model contains a free parameter R which determines the curvature of physical space and, as we shall see, has 
crucial effect on the number of solutions. This is a toy model which mimics critical phenomena in formation 
of solutions in the Einstein-Skyrme model 0 (but not only) and in principle it was the main motivation of 
the present work. 

The hedgehog Ansatz leads to solutions which in the original Skyrme model were called skyrmions. In 
general we call skyrmions all classical static solutions of the model which minimize energy within a given 
sector of the baryon number B. Equations of motion for the Skyrme fields are the Euler-Lagrange equations 
for the energy functional and, as such, the solutions may be not minima but saddle points. 

In the case of the Skyrme model on flat space, if U tends to a certain limiting matrix at spatial infinity 
independent of any direction, we can define U{ oo) = limiji^oo U(x). If so, we can attach to R 3 the point at 
infinity and by that complete the space to a compact one. Next we can identify points from R 3 U {oo} with 
points from S 3 by a suitable mapping - for instance the stereographic projection. A mapping from R 3 U {oo} 
to S 3 and further from Tx (R 3 U{oo}) to Tx5 3 is not an isometry and produces space-time with metric tensor 
whose curvature is not vanishing. It means that the Skyrme model on flat space is not metrically equivalent 
to the transformed one but topologically is. By means of that topological equivalence we can examine certain 
topological aspects of the flat Skyrme model in compact space such as S 3 what partly explains my work on 
the Skyrme model on S 3 . But we can look at this from quite another point of view. Physically a space-time 
endowed with non flat metrical tensor stands for gravity but the space-time from this point of view is then 
only a scene for skyrmions; space-time influences their motion but they do not influence its geometry. It may 
of course be so only under certain conditions e.g. if the presence of skyrmions only slightly disturbs metric 
tensor (like gravitational waves far away from their source). Then equations of motion for skyrmions decouple 
from Einstein’s equations and may be solved separately. 

My aim here is not investigation of the full set of equations of the Einstein-Skyrme model (see 0). I am 
looking for some aspects of nonlinear field theory coupled to gravity independent on the specific items of 
the space-time. The Skyrme model on S 3 is a toy model and as such aspires only to give insight into these 
unexplored areas of nonlinearity under gravitation. Here we treat the model only as an example of nonlinear 
field theory. 


2. The hedgehog Skyrme model on S 3 


In this paper I consider only static solutions of the Skyrme model on three-sphere which are critical points 
of the energy functional (1.2). Using the hedgehog Ansatz (]1.4|) we get 


E[F\= J 47T sin 2 ■i/idV’ 
o 


Rf 


2 r 


F' 2 + 2 


sin 2 F 
sin 2 ip 


+ 


sm 


2 F 


Re 2 sin 2 ip 


F' 2 + 


1 sin 2 F 

2 sin 2 ip 


In the limit e —> oo we regain the energy functional (with the unit of energy Rf 2 p/2) for harmonic maps between 
three-spheres which are equivariant with respect to the action of SO( 3) group and which were examined by 
Bizon in 0. Thus the Skyrme term may be considered as a perturbation which deforms harmonic maps of 
Bizon. 

Two parameters e and f n (in units in which c = 1) provide us with the natural unit of length L a := (e/^) , 

determining characteristic size of a soliton, and with the unit of energy E a = f 7r e~ 1 / 2. The second scale of 
length L g is provided by the radius R of the physical space - the three sphere. Roughly speaking, R artificially 
introduces gravity. (In the Einstein-Skyrme model, the second scale is introduced by the combination sfG/e 
originating from the full Einstein-Skyrme action §)• These two scales of length constitute a dimensionless 
parameter L := L g /L a = ef^R which is a free parameter of the Skyrme model on three-sphere. When e and 
f„ are nonzero and fixed, then L has the interpretation of a dimensionless radius of the three-sphere. In these 
units the energy functional takes the form 


TV 

(2.1) E[F] = 47 tL J sin 2 ip dtp 


F 


i/2 


^sin 2 F 
sin 2 ip 


2 

L 2 


F' 2 + 


1 sin 2 F 

2 sin 2 ip 


sin 2 F 
sin 2 ip 


Critical points of the functional ©> are solutions of the Euler-Lagrange equation 
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This is the quasilinear second order equation for the shape function F(ip). In what follows we shall refer to 
the parameter k 2 = 2/L 2 as the coupling constant. The only nontrivial solution of this equation, which is 
known analytically, is the identity map F(ip) = ip which has energy 

(2.3) E = 67r 2 (L + L _1 ) > 12 tt 2 ; 


where, depending on whether L > 1 or L < 1, the harmonic or the Skyrme term prevails. 

We are interested only in solutions for which both F(0) and F(ir) are integer multiples of n. Another 
possibility are the solutions for which F( 0) or F(ir) (or both) are some odd multiples of 7 t/2 (see the appendix 
for appropriate proofs). These solutions behave like the function \/Asinln(A) about A = 0, oscillating infinitely 
many times in the vicinity of poles. From the physical reasons we give up considering them, since otherwise 
the energy integral (HD would be divergent. 

For solutions, for which both F( 0) and F(n) are integer multiples of 7r, the Bogomolnyi bound holds, and 
we get it from (^d]) by completion into the squares 


7T 

(2.4) E[F] = AnL J sin 2 ^di/>| 


pi _ 1 sin 2 F 
L sin 2 if) 


n 2 2 sin 2 F 

+ 


L 2 sin 2 i/> _ 


f(tt) 

(F'-L) 2 + 3LF' 1 > 24tt J sin 2 FdF = 12n 2 Q. 


This bound is quite general and holds for all SU( 2) Helds of the model ( |l.2| ). It is clear that the Bogomolnyi 
bound is saturated by the identity map at L = 1. By the way the solution which is the counterpart of the 
flat one-skyrmion, converges uniformly to the identity when L \ y/2. The bound 12-7T 2 is not attained in the 
case of flat skyrmions. 

By equatorial reflection symmetry of the energy functional ( |2.l| ), which does not change energy, we can 
assume, without loss of generality, that for solutions F(n) — F( 0) ^ 0. For solutions the integer 


Q = 1(F(7t)-F(0)) 

is the topological charge which is topological constant of motion and cannot be changed by continuous 
deformation of a solution. 

Because of symmetries of the functional m we can narrow down the group of solutions. It suffices to 
take into account only solutions for which F( 0) = 0 and F'( 0) ^ 0, for if F(ip) with F( 0) = kn, 4 e Z is a 
solution then F{ip) — kn and —F{ip) are solutions too. Moreover if F(ip) is a solution with topological charge 
Q and with F(0) = 0, then G{ip) := Qn — F{n — ip) is a solution with G(0) = 0 and with the same topological 
charge Q but with G'(ip) = F'(ir — i/>). Hence it suffices to consider only solutions for which F'(0) A F'(n). 

So, to sum up, we can limit ourself, not losing generality, taking into account only solutions for which 
F( 0) = 0, F'(0) ^ 0 and F'( 0) ^ F'{n)\ of course in special cases F\n) can be negative. Numerics shows 
that for such solutions Q ^ 0. 


3. Tripartite structure of solutions 


3.1. Numerics. The aim of this work is to find the whole spectrum of solutions of the nonlinear equation 


(2.2). This can be achieved only by finding numerical solutions. For this purpose it is essential to know the 
Taylor series expansion of solutions in the vicinity of singular points of the equation, namely about ip = 0 
and ip = tv. The formal Taylor series expansion for solutions with F( 0) = 0 is 


(3.1) 


F{ip) = aip — 


i(a 2 — 1)(4 + a 2 n 2 ) 


ip 3 + 0(ip 5 ), a = F'{ 0). 


30(1 + a 2 n 2 ) 

This series expansion can be used as a launcher from the singular point to the domain of regular points 
of equation (2.2) where numerical integration may be successfully carried on. In the left neighbourhood of 
ip = 7r the series is similar. By the method of ’shooting to a fitting point’ Q we are able to get effectively 
numerical solutions of the equation. The series quite accurately approximates solutions at the boundaries 
and the n + 1 term can be used to assess an error one does by terminating the series at n’th term. This allows 
us to move away as far as possible from the singular points, within the accuracy desired, and then carry on 
the integration numerically. 

We denote by F 0 the function which we get by carrying on an integration forward step by step via equation 


(2.2), propagating the initial data F( 0) = 0 and F'(0) = a 0 from the north pole of the base three-sphere. 
Similarly by F 7r we denote the function derived by propagating the initial data F(n) = Qn and F'(n) = a 7r 
backward from the south pole. By assumption a 0 is positive. The numbers a 0 and we call the shooting 
parameters. 

Within a given topological sector with topological charge Q ^ 0 and for fixed nonzero radius L, there exist 
some finite number of pairs {a 0 ,a,r} of critical shooting parameters, such that the functions F 0 and F nr and 
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their first derivatives, launched from the opposite poles, will meet at some fitting point ipf which, from the 
symmetry, is taken to be the equator ipf = In this way, for these exceptional pairs {a Q , a n }, we get two 
complement ary p arts of the same integral curve which can be matched together to form the global solution 
of equation (T2) with the initial data (0, a 0 ) at ^ = 0 or equivalently (Q 7r, a^) at ip = n. In special cases 

d-jY — do Dr d-]{ — do 0. 

In generic situation if we start numerical integration with F(0) = 0 and with an arbitrary shooting 
parameter a 0 , then there will exist an integer l such that the numerical solution F„ pip ), continued up to 
ip = 7r, will attain the value F 0 (n) = (5 + 0 71 "- I n this way we get the function y[d 0 ) = F 0 (n) which is 
antisymmetric with respect to a Q and can be restricted to the set d 0 £ M. + . The function takes values in the 
set X := {1 £ K.: x = (i + 1)tt, l £ Z}, apart from some discrete number of critical values of a 0 for which the 
function attains the values being some integer multiples of n. Since y is a step function, if a 0 is critical then 
the difference y(a 0 + rf) — y{a 0 ) = ±7t/2, for sufficiently small y, is positive or negative in dependence on the 
value of d 0 , and for fixed a 0 on the sign of y. Having found numerically the function y(d 0 ) at some fixed radius 
L , we are able to find all critical shooting parameters a D . Of course this can be done only approximately 
since, numerically, one can generate the function y[d 0 ) only for a discrete number of values of a 0 ■ It should be 
clear that among the (approximately) critical values there may exist pairs which are (approximate) shooting 
parameters of the same global solution. In particular, if a global solution has a reflection symmetry in the 
base space, its pair of critical shooting parameters may consist of {a 0 ,a 0 } or {a 0 , — a D }. Next the pairs of 
approximate shooting parameters can be used as trial ones to find exact values for {a 0 ,a T } and afterwards to 
find the global numerical solutions of equation (2.2). It is achieved by zeroing a norm of the vector function 
( F a — F„, F' a — Fp) of two variables a 0 and a w , taken at the fitting point ipf = 7r/2. For this purpose it was 
used the globally convergent Newton-Raphson method for nonlinear systems of equations M. 


3.2. Terminology and notation. Numerics shows that all solutions of the Skyrme model on S 3 in general 
consist of three clearly distinguishable parts: a harmonic map in the middle and two skyrmions, each with 
certain topological charge, attached at the ends of the map. In other words in the vicinity of poles of the 
base three-sphere the solutions manifest their skyrmionic nature and in between they contain a harmonic 
map known from the paper of Bizon [fil. This should be understood in the following way. 


Firstly we consider skyrmionic ’arms’ that is the integral curves of equation (2.2) in the vicinity of poles. 
We focus our attention on the north pole. (In the case of the south pole one proceeds in the analogous way). 
In order to pass from the Skyrme model on S 3 to t he ordinary Skyrme model on flat space, we extract the 
skyrmionic part from the equation of motion (T2) introducing new dimensionless variable r := 2 Lip and 
defining the function Spr) := Fpip (r)). In a small neighborhood of the north pole r -C 2 L and consequently 
sin ip w r/2L. For all finite r we formally get, in the limit of infinite radius L (the unit of length is fixed), the 
equation for the limiting function Spr) 


(3.2) 


(\r 2 + 2 sin 2 S)S" + \rS' + S ' 2 sin 2S - \ sin 2 S - 


Sill 


2 S 


sin 2S = 0 


which is the equation for shape function of the hedgehog skyrmions on flat three-space jnj. Here r is the 
radial coordinate and the unit of length is e~ 1 Fp 1 with F n = 2f 7r . The difference in units comes from 


the conventions assumed by different authors. In the model, described by equation (3.2), for the reasons 
similar to those above, one looks for solutions which interpolate between 0 at r = 0 and nir at r = 00. Such 
solutions have finite energy and asymptotically behave as S n (r) = a n r + Opr 3 ) in the vicinity of r = 0 and 
as S n (r) = mr + (3 n /r 2 +o(r~ 2 ) at infinity ||]. Comparing the first series expansion with (HI) it is clear that 
^'(0) for skyrmionic arms must behave like 1/L in the limit L —> 00. We will denote by s„ such limiting arm, 
shrunk into the north pole, where n is a topological charge which simultaneously is the topological charge we 
ascribe to the corresponding skyrmionic arm of a solution of equation (2.2) when L is finite. The skyrmionic 
arm of the last solution we will denote by S„ to indicate that it evolves into s„ as L —> 00. For the south pole 
we proceed in the similar way putting r = 2LpK — '!/’)• I n both cases (the north and the south pole) we have 
to do with skyrmions if F'( 0) or F' pit) diverge to +cxd or with anti-skyrmions if the respective derivatives 
diverge to —00 as L —> 00. 


In the second step we consider another case in which for a solution of equation (2.2) in the limit L —> 00 
F'{ 0), F'(tt) or both are finite. If in this limit F'( 0) (or F'(n)) is finite then the respective skyrmionic arm is 
absent and we have to do with a harmonic map directly attached to the north pole (the south pole) or with 
pure harmonic map if both derivatives are finite. Harmonic maps appear in our model since if one drops the 
Skyrme term, what formally can be done by taking the limit L —> 00 in equation ( |2 . 2Q , we regain the equation 
for equivariant harmonic maps between three-spheres. These maps were first found and their existence proved 
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by Bizon in |lj. They are solutions of the equation 

(3.3) F"(ip) sin 2 ip + F' (ip) sin 2ip — sin 2F(ip) = 0. 

It was proved in [|lj that smooth maps F : S 3 —> S 3 , which are solutions of above equation, contain two 
countable families of harmonic representatives in the homotopy classes of degree 0 and 1. For a given map 
the index k of a map is defined as the number of times the map crosses the line F = tt/2. The topological 
charge is unit for odd and null for even maps. We will denote further these maps by the symbol h k . Solutions 
of equation ( |2.2| ) which in the limit L —> oo tend uniformly to h k will be denoted by H k and, for convenience, 
called harmonic maps, though in fact are not harmonic maps. 

To sum up saying harmonic or skyrmionic part of a solution of equation (2.2) we mean the limiting 
behaviour of the solution as L —> oo and use these names for convenience and for the sake of clarity of further 
considerations even if L is finite, though for finite L there is no explicit boundary between harmonic and 
skyrmionic part. Put differently in different regions of physical space - the three-sphere - solutions manifest 
their different nature. From now on by saying ’a solution is composed of’ a harmonic map or a skyrmion, we 
will mean that the solution in the limit of large L solves approximately - between the poles - the equation 
of harmonic maps ([hi]) and that in the vicinity of poles (after appropriate rescaling of variables as above) 
solves approximately the equation known from the flat Skyrme model ( [>. 2| j . We shall use above terminology 
regardless of a value of the radius L. 

Because of symmetries of the model we can identify solutions which can be transformed one into another 
using translations and reflections in the target space and reflections in the base space. From now on we will 
call such class of equivalence a solution. Each class has its own representative for which F'(0) > 0 and Q ^ 0 
simultaneously, so that by saying a solution we mean its representative. Each solution of ( |2.2| ) we denote 
by the symbol S n H k S m which informs us that in the limit of infinite radius of the base space, the solution 
takes the following limiting form s n h k s m : skyrmions s n and s m with topological charge Q = n or Q = m 
localized respectively at the north or the south pole of the base space, and a harmonic map h k with some 
index k which connects them. We admit of k = 0 and then H 0 denotes the neutral solution which is the 
representative of trivial solutions F = krr. Depending on whether the harmonic map has an odd or an even 
index, it contributes or not the unit topological charge Q h = sgn(F'(0)), where F is a solution of equation 

(3.3) . Thus the limiting solution s n h k s m has the topological charge Q = n + Qh + rn which of course must 
be the topological charge of the solution S„/-/ fc S m , since topological charge cannot be changed by continuous 
deformation within a class of functions fulfilling given boundary conditions. By definition for all solutions 
of equation (2.2) n ^ 0 (since F'( 0) > 0) and Q ^ 0. In special cases k = m = 0 and, in the limit of 
infinite radius, we have a pure skyrmion (localized at the north pole), or n = m = 0 and then we have a 
pure harmonic map h k = lim^—^ H k between three-spheres. A skyrmionic solution with negative topological 
charge we call anti-skyrmion and denote by S i.e. if m < 0 then S m = 5_ m . By h k we denote the harmonic 
maps for which F'(0) < 0. 

Profile functions of solutions of equation (2/2) and their phase diagrams are shown in this paper using the 
conformal variable x which is related to the angle ip by the equation x = In tan (^). The profile functions are 
plotted using coordinates (x,F(x)) and the phase diagrams using coordinates (F(x), F'(x)). 

Figure |l] illustrates above terminology using as an example the solution S 1 H 3 S 2 together with its phase 
diagram for different sizes of the base space . This solution, whose representant is S 2 H 3 S i, bifurcates from 
the solution S 1 H i S 1 at the critical value L « 126.56 of the radius of the base three-sphere. The greater 
the radius the more skyrmionic branches of the solution flatten and become tangent to the lines F = it and 
F = 2ty respectively for 1-skyrmion in the vicinity of the north pole and for anti 2-skyrmion in the vicinity 
of the south pole. The inner part of the solution in the limit of infinite radius overlaps the harmonic map h 3 
for which F(0) = tt and F(n) = 2ir. 

In the following we shall show that the extremely rich structure of solutions of the nonlinear equation (fh^), 
undergoes a strikingly simple scheme which enables us doing predictions of the properties of the solutions 
and how do they appear. 
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FIGURE 1 . Solution S1H3S2 and its phase diagram for different sizes of the base space. This 
solution whose representant is S 2 H 3 S 1 bifurcates from SiH 4 Si at some critical value of the radius 
of the base three-sphere. 


4. HOW DO THE SOLUTIONS APPEAR? 


In this section we assume that the parameters e and f w of the model are finite, nonzero and fixed once 
and for all. Using this parameters we define the units of energy f 7r e _1 /2 and of length f~ l e~ l . Then 
L = f^eR has the interpretation of a dimensionless radius of the base three-sphere. I emphasize the trivial 
fact the ’coupling constant’ n 2 = 2 /L 2 of the model can be changed both by changing physical radius R or 
by changing e or f n . In particular the limit k 2 = 0 can be attained either by taking e —> 00 , then we recover 
the sigma model with the unit of energy Rf^/2 (energies of harmonic maps are then finite and the same as in 
[|lj ); or by flattening the physical space putting R —> 00 , then we recover the flat Skyrme model with the unit 
of energy / r Tr e _1 /2 (and then the energies of skyrmionic solutions coincide with those from the flat Skyrme 
model). As we shall see this fact partly explains why the Skyrme model on three-sphere admits existence 
of branches of solutions of harmonic, skyrmionic or mixed type, which meet at some critical value of n 2 and 
disappear. 


4.1. Solutions which exist for all values of the coupling constant. For all n 2 there exist solutions 
which in the limit k 2 —> 0 tend to configurations H 0 , S 1 H 0 S 1 , S 2 H 0 S 2 , ..., S n H 0 S n , ..., with topological 
charge Q = 2n and solutions which in this limit tend to configurations H 1 , S 1 H 1 S 1 , S 2 H 1 S 2 , ..., S n H 1 S n , ..., 
with topological charge Q = 2n + 1. They all are invariant under the transformation F(x) —> Qn — F(—x) 
which for integer Q is a symmetry of the equation of motion (and in general is broken by generic solutions). 
As is seen on phase diagrams (fig.||) for small L there is no qualitative difference between solutions S^HiSn 
and S n H 0 S n except for the number of maxima for the function F '. However, in the limit of infinite L, they 
are quite different - S n H 1 S n contain additionally the harmonic map Hj^. This is reflected especially in the 
asymptotic behaviour of energies of solutions with odd and even topological charge in the two limits (figji]). 
In the limit f v —> 0 (no sigma term) Fty) solves the equation 


(4.1) 


F + F cot F - 


1 sin 2 F 

2 sin 2 if) 


= 0 


which can be formally obtained by taking the limit k 2 —> 00 in equation (2.2). As is seen on phase diagrams 
(fig#, in the limit L —> 0 the maxima of F'(i/)), localized at such ipk for which F(ipk) = kn where k is some 
positive integer, tend to infinity. This follows from the formal Taylor series expansion of solutions of equation 
( 0 ) whose singular part for F' in the neighbourhood of t/jk is 


In the limit L —> 0 the shooting parameters for the solutions tend to finite values e.g. 2.5635, 4.1047, 5.641, 
7.1755, 8.7089 respectively for S 1 H 1 S 1 , S 2 H 1 S 2 , S 3 H 1 S 37 S 4 H 1 S 4 , S 5 H 1 S 5 . These values form a sequence 
which asymptotically is arithmetic. The same (but with different values) holds for solutions S n H 0 S n . For L 
small enough the energies Eq of solutions S„H 1 S„ and S„/-/ 0 S„ behave like eg/L where constants eg form 
an ascending sequence in raising topological charge Q and asymptotically scale with Q 2 that is 


lim LEq(L) = eg, and lim —tt = const. > 0. 
L —>0 w < 5^00 Q 2 


These energies would become finite again if, instead of f n e 1 /2 as the unit of energy, we chose R x e 2 for R 
and e fixed (but then we would get quite different model). In the limit of large L the asymptotic behaviour 
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FIGURE 2. Solutions SiH 1 Si ( k, 2 = 10 —5 ,..., 10 3 ) and S 2 /-/ 0 S 2 (k 2 = 10 —3 ,..., 10 2 ) and their 
phase diagrams for different sizes of the base space. The smaller the radius L the higher the maxima 
for F'. In the limit L —> 0 the maxima become infinite and then the solutions tend to the solutions 
of the limiting equation (0). 



FIGURE 3. Energies (in units f^e^ 1 / 2) for solutions S„/-/ 0 S„ or S„HiS„ for n ^ 1. Counting from 
the bottom the curves are H 1: Si/-/ 0 Si, S 1 H 1 S 1 , S 2 S 2 and so on with growing topological charge Q 
up to S 4 H 1 S 4 . It should be evident how this curves will behave for Q > 9. The lowest line - the 
identity H 1 - can be explicitly calculated from equation (^^) hence we get the asymptotic behaviour 
of energies for other solutions. 


of energies for solutions S„H 0 S„ is very different from the behaviour of S^Hj^Sn. This difference comes from 
the fact (which is the numerical observation) that energies of n-skyrmions are bounded while the energies 
of harmonic maps are unbounded in the limit L —> 00 (with / r 7r e~ 1 /2 as the unit of energy). Due to the 
tripartite structure of solutions and the fact the integral of energy is additive, this behaviour is clear. This 
fact (clarified on figJJ) is quite general, that is if we take a solution S n H k S m then 
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FIGURE 4. This figure shows how the hypothesis of tripartite structure of solutions works in 
practice. The lower line shows the energy (multiplied by 2) of 1-skyrmion Si which, having bifurcated 
from the identity solution Hx at L = \/2, evolves becoming true 1-skyrmion (with energy « 1.23145- 
12-7T 2 ) localized at one of poles in the limit of infinite radius of physical space - the three-sphere. The 
other lines are: the energy for the solution Si/-/ 0 Si composed of two 1-skyrmions, and the energy 
of the solution Si Hi Si (which contains harmonic map Hi inside and 1-skyrmions attached to its 
endings) after subtracting the energy of the harmonic map Hi. 


and 


lim 


E Sl 


H k S, 


XL) 


= 1 . 


E Hk (L) 

From the tripartite structure it is also clear, that in the limit L —> oo, the s hoo ting parameters for solutions 
S„H 0 S„ and S rl H 1 S rl divided by L are the same as for n-skyrmions (see pt. 4.3). 

We can assume that the arms of solutions S„H 0 S„ and S n H 1 S n pass from the regime described by equation 
(0 to the regime when they are skyrmionic arms (the n-skyrmions), if their phase diagrams flatten at x = 0. 


4 . 2 . Harmonic maps and accompanying solutions. There exist a descending sequence of critical values 
of coupling constants fc ^ 2 at which there appear harmonic maps with index k. The only exception 

is the identity map H 1 which exists for all values of n 2 . In the limit k 2 —> 0 all these maps tend uniformly 
to those found by Bizon in [Q. These limiting maps are shown in figure |j. In the case when this limit is 
understood in the sense e —> oo (with Rf\f 2 as the unit of energy) the energies of limiting harmonic maps 
coincide with those from |l[]. If the limit k 2 —> 0 was understood in the sense R —> oo (with f„ and e fixed), 
then these energies (with / r 7 r e ” 1 /2 as the unit of energy) would be infinite. The energies of harmonic maps 
and the accompanying solutions are shown in fig.||. The energy of the identity solution H 1 can be computed 
explicitly from ( {2.l| ) 

E Hl (L) = 6 t t 2 {L + L" 1 ) > 12 t r 2 . 

This energy is infinite in the limit L —> oo and attains its minim um 127 t 2 at L = 1 which is the global 
minimum of energies as follows from the Bogomolnyi bound ( [2.4|) . Numerics shows that in this limit the 
energies of harmonic maps scale with L and that the functions En k (L)/(12 tt 2 L) have finite limits 


lim 
L —>00 


E Hk (L) 

12tt 2 L 



where the sequence {Ek}ke N (in units in which Rf 2 K /2 = 1) is ascending and converges to 2/3. This limiting 
value is the energy of the singular map F = 7r/2 to which the sequence of functions converges pointwise 

(but not uniformly because the singular map does not fulfill the boundary conditions). In the limit e —> oo 
the respective sequence of energies Ek is precisely the same as in the sigma model on S 3 examined in (TJ. 
Furthermore, as we shall see below, after the above rescaling E —>■ e = E /X, the energies of skyrmionic 
solutions tend to zero so in the limit of large L different solutions, with different topological charges, but 
containing the same harmonic map inside, have the same limiting ’energies’ e. Moreover their profile functions 
tend pointwise but not uniformly to those of harmonic maps they contain. This observation qualitatively 
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FIGURE 5. Even and odd harmonic maps H k for k = 2,4,6,8,10,12 and for k = 1,3,5,7,9,11 
respectively and phase diagrams for k = 12 and k = 11 in the limit of infinite e. In the limit k —> oo 
phase diagrams have infinite self-similar structure. 



FIGURE 6. Energies (in units where Rf\/ 2 = 1) for the first five harmonic maps and the associated 
solutions. Solid lines shows harmonic branches H k (only five for clarity), dotted lines - branches of 
associated solutions SiH fc _ 2 Si and SiH fc _ 2 Si for odd and even k ^ 2, respectively. These branches 
coalesce in a cusp for descending L. The short dotted lines are branches of solutions SiH k -i bi¬ 
furcating from harmonic maps H k . The identity map H i is exceptional and exists for all L. The 
solution Si bifurcates from Hi at L = y/2. 


explains the conjecture made by Y. Brihaye and C. Gabriel in ||, who found several solutions, that (using 
my terminology) 


K —* K 2k 

lim E, 


lim E„ 2 k {K) = K hm E Sl J7 2fc _ 2 s l( «) = E Hak (K 2k ), 


lim 


i™0®SiH 2fc Si( K ) = = E H2k (0), 


_« 2fc+1 ^ 2k+1 (") = „ Jlf k+1 E Sl H 2k _ lSl (*) = E H 2k+1 (« 2 fc+ l), lim E Si H 2k iSi (*) = lim E H2k _ x (*) = E^ (0). 

The numerical evidence for these limits, which are made here far more plausible, is shown in fig.^j and |i]. 

It should be clear, that in the limit of large L (understood as large e), the Skyrme term can be treated 
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FIGURE 7. This figure, which is asymptotically self-similar, shows branches of solutions for families 
of harmonic maps. The lines which asymptotically (k 2 —> 0) become horizontal are branches of 
tangents F'( 0) for the harmonic maps H k at ip = 0. Tangents F'(n) at ip = n are equal to F'( 0) and 
—F'( 0) respectively for odd and even harmonic maps. Asymptotic values of F'( 0) are precisely the 
same as for harmonic maps. If k 2 increases then each harmonic branch H k turns back at a critical 
value k\ and becomes a branch of the associated solution SiH fc _ 2 Si for k even and SiH fc _ 2 Si for 
k odd. All associated harmonic branches, in the limit k 2 —► 0, tend to the common branch being 
the line of tangents F'( 0) of the 1-skyrmion. The only exception is the identity Hi which exists 
for all k . It is the first mechanism of how the solutions are born. To the second mechanism there 
depend branches crossing the harmonic ones. The upper branch of the second mechanism is the 
line of tangents F'( 0) which asymptotically attains the branch of the 1-skyrmion, while the lower 
branch is the line of tangents F'(n) or —F'(n) if it crosses an odd or an even harmonic branch. 
Asymptotically the lower branch is tangent to the harmonic branch H k ~ i, so the solutions of the 
second mechanism have to be nonsymmetric solutions SiH k ~i■ If k = 1 it is 1-skyrmion alone whose 
lower branch is asymptotically tangent to the harmonic branch H 0 - the vacuum solution F = 0. 


as a perturbation which deforms harmonic maps from the sigma model on S 3 . On the other hand the limit 
of large L has its unperturbative regime which is attained by fixing e and taking the limit R —> oo, thus 
there appear other solutions accompanying harmonic maps. There exist two types of such solutions which 
cannot be derived by perturbation from harmonic maps. A solution of the first type appears together with 
accompanying harmonic map in a characteristic cusp on the diagram of energies at critical radius L. Next 
they evolve on separate branches. (By a branch we mean the plot of the energy or of the shooting parameters 
of a solution as a function of L or k 2 ). A solution of the second type formes another branch which bifurcates 
from the harmonic map’s branch elsewhere at a larger critical radius L, but before the next harmonic map 
with higher index appears (see fig.|7j). These two mechanisms of appearing of solutions, namely the creation 
of a pair and next the bifurcation of a second solution from the just created harmonic branch, as we shall see 
further, is generic and will be repeated during appearing of other solutions with higher topological charges. 
The mechanisms in the case of families of harmonic maps are described below. 

Mechanism I: Each harmonic map H k , k = 2 , 3 ... appears together with an accompanying solution. 
The pattern of appearance is as follows. The even harmonic map H 2k appears together with the 
solution S 1 H 2k _ 2 S 1 e.g. H 2 together with Sib/oSj while the odd harmonic map H 2k+1 appears together 
with S 1 H 2k _ 1 S 1 e.g. H 3 together with S 1 H 1 S 1 . See tab. [I] and fig. || 

The higher the index of a map, the larger L at which it appears. Here again, in connection with harmonic 
maps, as in [l]] there appears the characteristic number e 27r /^ but this time it is the limit of a sequence 
defined in the following numerical hypothesis: 

Hypothesis 4 . 1 . The sequence of quotients of critical coupling constants at which there appear 

odd or even harmonic maps is decreasing and attains the limiting value (~ 10.749087) 


lim = e 27r/x/ ^ 
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mechanism 1 

mechanism II 

pair 


Kfc/ Kk +2 

pair 


Kfc/Kfc+2 

h 2 

SiSi 

2.613327- 

10" 1 

— 

Si Hi 

HiSi 

5.245019 

• 10" 2 

— 

h 3 

SiHiSi 

1.706483- 

10" 2 

- 

SiH 2 

H 2 Si 

5.703710 

• IO " 3 

- 

Hi 

SiH.Si 

1.593848- 

IO " 3 

12.80482 

SiH 3 

H 3 Si 

4.372362 

• 10 -4 

10.95256 

H 5 

SiH.Si 

1.391083- 

10 -4 

11.07578 

Si Hi 

H 4 Si 

4.849196 

• 10" 5 

10.84536 

H 6 

SiHiSi 

1.351911- 

IO " 5 

10.85799 

SiH 5 

H 5 Si 

3.767585 

• 10 -6 

10.77275 

H 7 

SiH 5 Si 

1.197673- 

10 -6 

10.77724 

SiH 6 

H 6 Si 

4.188957 

•IO " 7 

10.75925 

H 8 

SiH 6 Si 

1.167865- 

io - 7 

10.75914 

SiH 7 

H 7 Si 

3.259313 

• IO " 8 

10.75149 

h 9 

SiH 7 Si 

1.036053- 

IO " 8 

10.75172 

SiH 8 

H s Si 

3.624800 

• IO " 9 

10.75006 

H 10 

SiKSi 

1.010587- 

10 -9 

10.75002 

SiH a 

H a Si 

2.82075 • 

10- 10 

10.74931 

Hu 

SiH a Si 

8.966418- 

10- 11 

10.74933 

Si Hi 0 

HioSi 




h 12 

SiH 10 Si 

8.746281 • 

10- 12 

10.74918 

Si Hu 

HuSi 

too poor 

numerical accuracy 

h 13 

SiHuSi 

7.760219- 

10- 13 

10.74911 

SiH 12 

H 12 Si 





Table 1. Solutions accompanying harmonic maps appear in two mechanisms. In the first 
mechanism there appear even harmonic maps H 2k together with S 1 H 2k _ 2 S 1 or odd harmonic 
maps H 2k+ 1 together with S 1 H 2k _ 1 S 1 . In the second mechanism new solutions bifurcate from 
already existing harmonic maps i.e. S 1 H 2k _ 1 , simultaneously with H 2k ^ 1 S ll bifurcate from 
H 2k while S t H 2k , simultaneously with H 2k S 1 , bifurcate from H 2k+1 . 



FIGURE 8. This figure illustrates the evolution of even and odd harmonic maps together with 
accompanying solutions. The profile functions on the left and respective phase diagrams on the right 
are plotted for different values of coupling constant. The evolution is exemplified by pairs of solutions: 

Hi together with SiH 2 Si at the top and H 5 together with SiH 3 Si at the bottom. Dotted lines shows 
the limiting solutions in the vicinity of critical values of coupling constant n 2 Hi = 1.593848 • 1CF 3 and 
Kh 5 = 1.391083 • 10 -4 . 

Using this hypothesis it is easy to predict when and together with what solution other harmonic maps will 
appear. This reveals the whole structure of appearance of harmonic maps together with their companions. 
But that is not the whole story. Numerics shows that there have to exist some (conformal) instabilities in 
solutions containing harmonic maps which lead to new solutions which bifurcate from the previous. This 
leads to the existence of the second mechanism of how solutions are born. 
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Mechanism II: from each even harmonic map H 2k there bifurcates the solution S 1 H 2k _ 1 and simulta¬ 
neously its reflection in the base space H 2k _ 1 S 1 . Similarly from each odd harmonic map H 2k+1 there 
bifurcates the solution S 1 H 2k together with H 2k S 1 which is derived from the first by reflection both 
in the base and in the target space. See tabj^ and fig jo]. 






FIGURE 9. The figure shows profiles and phase diagrams for solutions appearing in the second 
mechanism in which new solutions bifurcate from already existing harmonic maps. Upper figures 
shows Si Hi for several k 2 up to 10“ 6 (solid lines) and H 1 S 1 (dashed line); both bifurcate from H 2 
(dotted line) at k 2 = 5.245019-10~ 2 . The figures in the bottom show SiH 2 for several k 2 up to 10“' 

(solid lines) and H 2 S 1 (dashed line) - both bifurcate from H 3 at k 2 = 5.70371 • 10 -3 . It is easy to 
distinguish skyrmionic and harmonic components of solutions (e.g. see phase diagrams). In general 
from harmonic map H k k > 1 there bifurcates S 1 H k - 1 and simultaneously H k ~iSi for k even or 
H fc _iSi for k odd. In special case when k = 1 from harmonic map Hi there bifurcates 1-skyrmion 
Si at the north and at the south pole. 

Here again, as in the case of the first mechanism, one can form the sequence { of quotients of 

critical coupling constants at which new pair of solutions bifurcates from odd or even harmonic maps. Again 
this limit attains the magic value e 27r / %/ A 

It is worth noting that the index and topological charge of solutions appearing in pairs is the same. 
Similarly the same is in the case of an already existing solution and the one which bifurcates from it. I state 
this as an empirical fact which is true for all solutions of the Skyrme model on S 3 irrespective of a value of 
the topological charge, but by an index we mean then the index of a solution within the 7r-wide strip which 
contains the harmonic map. Summing up, the formation of solutions in the second mechanism is governed 
by the laws of conservation of the index and of the topological charge of appearing solution. 

4.3. The n-skyrmions and their companions. The solutions whose shooting parameters in the limit of 
large L have asymptotic behaviour of the type do ~ O(L) and simultaneously ~ 0(L~ 2 ), where a 0 and 
are the values of F'(0) and F'(tt) respectively, we refer to as n-skyrmions and denote them by S n where n 
is topological charge. The name skyrmion for these solutions is justified by the observation that in the limit 
of infinite radius L of the base three-sphere their energies are just the same as in ordinary flat Skyrme model 
(see tab]| and fig.[To|) . As we have already seen, after appropriate rescaling of independent variable, in the 
limit L —> oo the skyrmionic solutions fulfil the Skyrme equation ( |3.2| ) in flat space. These solutions we order 
with increasing topological charge which simultaneously is the order in which they appear as the radius L 
increases. A solution S„, for which n )> 2, appears together with its companion S n _ 1 H 1 which is composed 
of the harmonic map H, and attached to it (n — l)-skyrmion (see fig.[ll]). In the special case when n = 1 the 
scheme of appearance is quite different - the solution S 3 bifurcates from the already existing H, at the critical 
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Sr 

s 2 

S 3 

s 4 

s 5 

s 6 

s 7 

Ss 

^lim E n /(12n 2 ) 

^lim E n /(127v 2 n 2 ) 

1.23145 

1.23145 

3.66707 

0.91677 

7.33917 

0.81546 

12.2539 

0.76587 

18.4135 

0.73654 

25.8186 

0.71718 

34.4701 

0.70347 

44.3680 

0.69325 

lim ao tn /L 

lim ao, n /(ivn) 

2.00754 

2.00754 

3.93062 

1.96531 

5.77722 

1.92574 

7.59385 

1.89846 

9.39525 

1.87905 

11.1877 

1.86462 

12.9745 

1.85349 

14.7572 

1.84466 

lim L 2 a n n 

L—kx> 

lim L 2 a n n /n 2 

L—>00 ’ ' 

0.45830 

0.45830 

1.3707 

0.3427 

2.7540 

0.3060 

4.6096 

0.2881 

6.9386 

0.2775 

9.7350 

0.2704 

13.0132 

0.26557 

16.759 

0.2619 


Table 2. Limiting values of energies (with f 7T e ~ 1 /2 as the unit of energy) and scaling of 
shooting parameters of skyrmionic solutions. These limiting values scale again with topolog¬ 
ical charge in the limit of large values of topological charge. 



FIGURE 10. The energies (in units f- n e~ 1 /2) of solutions S„ and S„_iA/i for n= 1,2, ...8. The 
asymptotically horizontal lines are energies of skyrmionic solutions S„ which, in the limit of infinite 
radius L, agrees with the energies of skyrmions from the Skyrme model on flat space. The other 
lines, which join the branches of S„ in cusps, are the energies of the solutions S„_iA/i. The only 
exception is the solution H 4 which exist for all L. For L large enough the energy of H± behaves as 
6n 2 L. The same asymptotic behaviour holds for S„_iA7i. The larger topological charge Q = n, the 
larger the radius at which S„ appears. 

radius L\ = \/2. The table |j shows the sequence of critical radii {L n } n ^2 and respective accompanying 


pair 

Ln 

L n /n 

pair 

Ln 

Lnj Tl 

S 2 

S 1 H 1 

3.36368 

1.68184 

s 7 

SeHr 

12.0117 

1.71596 

S3 

S 2 H r 

5.1125 

1.70417 

Ss 

S 7 H 1 

13.7306 

1.71633 

s 4 

S3 Hr 

6.8442 

1.71105 

s 9 

Ss H 1 

15.4489 

1.71654 

s 5 

S.AL 

8.5694 

1.71388 

S 10 

SeAVr 

17.1666 

1.71666 

s 6 

s 5 h 1 

10.2914 

1.71523 

S 41 

S W H 1 

18.8839 

1.71672 


Table 3. Skyrmionic solutions S„ with n ^ 2 and their associative solutions S n - 1 H 1 appear 
at critical values of radius L n . The sequence { L n /n } n has a finite limit. 


solutions. 

Numerics shows that critical values L n scale with topological charge Q = n and supports the hypothesis 

that the limit of the sequence {L n /n } n6 n exists and is finite 

n . L n 
lim — = s c . 

n —kx) 77 , 

A fit gives the value s c ~ 1.7169. From the table [^. it is clear that shooting parameters ao and for 
n-skyrmions 5„ scale asymptotically with the radius L that is the limits limi—jocOo/L and limj^oo a v L 2 
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FIGURE 11. (a) The first four pairs of solutions S„+i and S„H i for n = 0,..., 3 at L = 100\/2 and 

(b) their phase diagrams. At this L there already exist skyrmions with topological charge up to 82! 

(c) The 1-skyrmion Si bifurcates from the harmonic map Hi at L = \/2 and evolves to its limiting 
configuration - 1-skyrmion from flat space; the evolution is shown for radii up to L = IOOOa/ 2 which 
form geometrical sequence. 


exist and are finite (see fig.O). In addition this asymptotic equalities, as well limiting energies, scale again 



2 

K 


FIGURE 12. Shooting parameters a 0 = F'{ 0) and a n = F'(n) for solutions S„ and S„_iHi, 
n = 1, 2,... 11. For the identity solution Hi a 0 = a n = 1 (the horizontal line) for all k 2 . For other 
solutions a 0 > 1 > a„. For all n > 2 the branches of a* and a 0 of solutions S„ and S„_iHi meet at 
critical value k^ and disappear, besides for all n < ^n-i an d for all k 2 < : o 0 (S„) > a 0 (S„_iHi) 

and a^(S„) < a^-(S„_iHi). Moreover, a 0 {S n -iH 1 ) tends asymptotically to a 0 (S„_i) and a^Sn-iHi) 
tends to a n (H 1 ) = 1. 


with topological charge and this all is summed up by the following limits: 


V 1 E n {L) 

lim lim - -r— = s e 

n— >oo L—>oo 1Z7T Tl 


y y a 0 ,n(-^) 

lim lim -= s 0 

n —>oo L—XX) Ln 


lim lim L 2 a7r ’"j —1 = s 


n—>• oo L—> oo 


Fits give approximate values to the above limits: s e ~ 0.688, s a « 1.77, ~ 0.24. The value of the 

asymptotic energy of the skyrmionic solution Si is the same as the energy of 1-skyrmion on flat space and 
confirms results given in || where the Skyrme model on S 3 was investigated. From the table ||. it is clear 
that n-skyrmions with Q = n ^ 2 are not energetically stable because E n /E\ > n. In particular, E-i > 2£’i 
what suggests that dynamically S 2 would decay into a system of two 1-skyrmions (maybe SiHoSi?). This 
suggest that n-skyrmions are unstable (for instance in the case of flat space skyrmions stable configurations 
are not spherically symmetric). 
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5. Full structure of solutions of the Skyrme model on S 3 

The main result of this work, besides the statement of the tripartite structure of solutions of the Skyrme 
model on S 3 , is the discovery in what configurations do the solutions appear. My analysis of a vast number of 
numerical solutions of equation ( |2.2[ ) and their evolution with the radius of the base three-sphere, had led to 
stinkingly simple picture. If a new solution appears when the radius exceeds its critical value, characteristic 
for that solution, then either it appears together with the second accompanying solution or it bifurcates from 
an existing solution. The clue in guesswork was to divide all the solutions into disjoint families n 3 ki as shown 
in the table [|. The nonnegative index n means that within the family n 3 ki we ask how the solution 5„ H p ( i k) S r , 


n 

n s k0 

iqk 1 


pairs of solutions 

Q 

pairs of solutions 

Q 


s„tf 2 fe _ 2 s„ s„_ 1 /-/ 2 fc s „_ 1 

0 

S n H 2k _ 1 S n _ 1 S„_ 1 /-/ 2fc S„_ 1 * 

0 

k> 0 

s„/-/ 2 fe _ 2 s„_ m s„_ 1 /-/ 2 fe _ 1 s„_ m 

m 

s„/-/ 2 fc _ 1 s„_ m s„_ 1 /-/ 2 fc s„_ m 

TO — 1 


s„/-/ 2fe 2 s „ 1 s„ 1 /-/ 2fe 1 s „ 1 a 

2n — 1 

S„H 2fc 1 S„ S„ 1 /-/ 2fc+1 S„ 1 

2 n — 1 

k = 0 

s„s„ 

2 n 

S„H 1 S n 

2,71 -f-1 


Table 4. Disjoint families of solutions n 3 k i of the Skyrme model on S 3 where n £ N U {0}, 
* = 0,1 (we say that i = 0 is opposite to i = 1). Solutions S„S n and S n H 1 S n which exist 
for all radii L of base three-sphere are included conventionally in the families. The solutions 
from njki marked by the star * belong to families with subscript n — 1 and with the opposite 
i and are unstable at certain critical value of L when they ’emit’ a solution written on its 
left side. The positive integer mn takes the values 1, 2,..., 2n — 2 if i = 0 or 2,3,..., 2n — 1 if 
i = 1. 


which contains the n-skyrmion at the north pole, harmonic map H with index p inside, and r-skyrmion at 
the south pole; appears under assumption that r takes such values that r does not exceed n and that the 
topological charge of the solution is not negative. The topological charge depends of course on whether the 
solution contains an odd or an even harmonic map. To see how this all works I give here several examples: 


Solutions S„H 0 S„ and SnHjS,,: The solutions S„H 0 S„ from mOO never decay in contrast to the 
solutions Snb/iSn from m01 from which S„ +1 S„ bifurcate belonging to the families n+d-,10. As 
we saw before the solutions S„S„ and S n H 1 S n exist for all radii of the base three-sphere. 
Solutions H fc (’harmonic maps’): Families l 3 k0 and l 3 kl contain even and odd harmonic maps 
respectively. From the table [|. we read off that the maps H 2 , A7 4 , ..., appear together with 
SiSi, S 1 H 2 S 1 , • ••, respectively, and odd maps H 3l A7 5 , ..., appear respectively together with 
S 1 H 1 S 1 , SxHsSi,.... We saw in §4.2 that the greater the index of a harmonic map, the greater 
the radius L at which it appears. Nevertheless, before H p+1 appears, S 1 H P _ 1 bifurcates from 
H p . Of course, the solutions H p for p = 2, 3,... are not harmonic maps (we have used this name 
for convenience), but in the limit of infinite radius of the base three-sphere they also fulfil the 
equation of harmonic maps between three-spheres, 
skyrmionic solutions: Families n 3 10 for n = 2,3,... contain n-skyrmions S„ which appear to¬ 
gether with S n _ 1 H 1 . We remember from how to interpret the name n-skyrmion of the solution 
S„ which, for finite radii of the base three sphere, is not n-skyrmion at all. But the similarities 
and proper interpretation introduced in previous paragraphs had justified the terminology. 
Solutions with Q = 0: contain even harmonic maps which can appear only in pairs S„H 2k S rl 
together with S n+1 H 2k _ 2 S n+ll whereas solutions containing odd harmonic maps S 7l H 2k _ 1 S„ +1 
appear from decay of S p H 2k S p where k ^ 1 and n ^ 0. Further it turns out that the greater 
n with fixed k or the greater k at fixed n, the greater critical radius at which given solution 
appears. 


Numerics shows that in general there is no rule which could be used to ’’forecast” the order of appearance 
of arbitrary solutions (at least if one uses numerical methods of finding solutions). If we take the family of 
solutions with n-skyrmion at the north and m-skyrmion at the south pole fixed but which are distinguished 
by the harmonic maps with different index k they contain, then the cusps of the type as in fig.[l(j| or knee-like 
shapes as in fig]?] and in fig.[l2], locate, in the limit of large k, along straight lines (if we use logarithmic 
scales) whose slopes are different in dependence on a pair (n, m). In this limit, along given straight line, 
these critical values are distributed homogenously what implies that asymptotically they form a geometrical 
sequence with base quotient being certain power of . Now it is clear that the sequences of critical 


















STRUCTURE OF SOLUTIONS OF THE SKYRME MODEL ON THREE-SPHERE 


17 


values from different lines are not commensurate with each other and that is why it is impossible to predict 
which solution would appear as the next one after the one arbitrarily chosen. 

In spite of above we can distinguish certain ’directions’ in which such ’putting in order of appearance’ is 
possible. Here are some of them: 

(1) If we take the set of solutions S„H p S m , with arbitrary but fixed n and to, with different p 0, then 
the larger p , the greater the radius at which the solution containing harmonic map H p appears. 

(2) For fixed topological charge, the solution S„ +s /-/ p S m+3 appears at larger radius than the solution 

S„H p S m . 

(3) For fixed nonnegative n, the solution S„ + 3 S m appears at larger radius than S„S m where s ^ 1 and 

0 ^ — to < 2 n. 

(4) If positive n is fixed then sequentially there appear S n S„_i, S„S„_ 2 , ..., S„5_„ (see fig-IIdl). 



FIGURE 13. Families of solutions from n,10 for n = 2,3,4. A number at each branch, showing 
shooting parameters a 0 as functions of k 2 , is the topological charge of skyrmionic solution localized 
at the south pole. 


We can include solutions containing harmonic maps into (js|) and (Q). Using above rules, reflections in the 
base and/or in the target space, and using the table we can generate much more rules. But they are not 
all of course. Using above empirical rules we are able to ’predict’ for example that the solution S 2 Si appears 
earlier (at lower L ) than S 2 ; it means that although there exists a solution ( e.g. S 2 SP there do not exist all 
its components yet ( i.e. S 2 does not exist yet). 
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Appendix A. 


A.l. Local existence theorem and asymptotics for solutions. In section (4.3) by change of variable 
ip = r/2L we defined the function Sl{t) = F (’ip)\ 1 f,— r / 2 L which in the limit L —> oo fulfils the Skyrme equation 
Solutions of the equation are critical points of the energy functional 


(A.l) 


E[S] = 4tt 


S" 2 + 2 


sin 2 S 


s'- 


1 sin 2 5\ sin 2 S 

2 r 2 


2 dr, 


where a 2 depends on a unit of length one chooses (the second integrand above is not scale invariant) and 
here a 2 = 8 . For proofs of asymptotic behaviour of solutions of the equation 

1 0 1 _ gi n 2 TP 

(A. 2 ) (sin 2 ip + k 2 sin 2 F)F" + sin 2 tpF' -\ — n 2 sin 2 F{F ') 2 — sin 2 F - n 2 — 5 — sin 2 F = 0, 

2 2 sin 1/1 
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the crucial point is the introd uctio n of a fun ction Q which (by analogy with can be simply guessed by 
comparison of the functionals (A.l) and dU- 

Definition A.l. 

(A.3) Q(xf) := -F ,2 (sin 2 xf + k 2 sin 2 F) - i sin 2 F - xf G [0 ,7r]. 


In this section I exploit methods elaborated by McLeod and Troy in their work || where they gave 
theorems concerning solutions of the Skyrrne model on R 3 , and in principle I repeat their ideas extending 
them on equation (A.2). In this (numerical) work I need only the proof that the series expansion ([hi]) for 
solutions of the Skyrrne model on <S 3 for which F(0) and F(tt) are integer multiples of 7r is valid. 

Lemma A. 2 . As xf \ 0 the function Q{xf) sin 2 xf converges to a finite limit. 

□ . Let xf G [0, 7 }Tt) then the function g = sin 2 xf(Q + L) is increasing because from the definition of Q 

g' = sin 2 xf(Q' + 2 Q cotxf + cot xf) = 4 sin 2 F cot t/>(|k 2 F' 2 — i) + cot xf > cos 2 F cotxf ^ 0. 

So that g decreases as xf \ 0 and the function Qsin 2 ^> has to attain a limit which is finite or not. But it 
cannot equal — oo what would imply from Q sin 2 xf that j sin 2 ^>F ,2 (sin 2 xf + k 2 sin 2 F ) was negative what is 
impossible. Hence there exists a constant a such that 


(A.4) 


lim Q(xf) sin 2 xf = a, and a is finite. 

■0—*-o+ 


□ . 

What is the consequence of this limiting behaviour of Q(xf) for convergence of F(xf) at xf = 0? The answer 
is the following lemma. 

Lemma A.3. As xf \ 0, the function F(xf) tends to a limit. 

□ . For contra dictio n sup pose that F(xf) does not tend to a limit as xf \ 0, hence sinF does not tend to a 
limit too. From (A3) and (A.4) it is seen that 


and next that 


\F' 2 sin 2 i/>(sin 2 xf + n 2 sin 2 F) — | sin 2 xf sin 2 F ~ a + g n 2 sin 4 F 
WF' 2 sin 2 xf sin 2 F ~ a + \n 2 sin 4 F 


where we used the fact that since sin 2 F does not tend to a limit then (sin 2 xf + n 2 sin 2 F) ~ k 2 sin 2 F and 
since sinF is bounded then — \ sin 2 xf sin 2 F —> 0. Moreover we can choose a descending sequence of points 
{xf„} in such a way that \\m. n ^ 00 xf n = 0 and that simultaneously fn 2 sin 4 F{xf n ) ^ —a (if a < 0) and 
i« 2 sin 4 F[xf n ) ^ 0 for all n. For these xf n jK 2 F 12 sin 2 xf sin 2 F is finite and different from 0. Having noticed 
that 

F' sin xf = xfF'{l +o(xf)) if xf < e, 

we see that xf n F'(xf n ) is finite and different from 0. In particular there exist such A that 0 < A < 1 and for 
all n xfF'{xf) is bounded and different from zero in some interval xf n ( 1 — A) ^ xf ^ xf n ( 1 + A). On the other 
hand 

( Q(xf) sin 2 xf)' = sin 2 F cot xf(^n 2 F' 2 sin 2 xf — sin 2 xf), 

and since xfF' is finite and bounded from zero there must exist such constant M > 0 that for sufficiently 
large n 

( Q(xf) sin 2 xf)' > A- if only xf n {\ - A) < xf < xf n ( 1 + A). 

tyn 

Integrating this in the interval xf~^xf^ xf+ where xf± =xf n { 1 ± A) we get that 

s (xfn) ~ s(xf~) ^ 2MA, and s(xf) = Q(xf) sin 2 xf. 

But this says that Q (xf) sin 2 xf as a function does not fulfils the Heine criterion for convergence, what contra¬ 
dicts the fact proved in previous lemma that Q{xf) sin 2 xf tends to a limit if xf \ 0. Thus F(xf) tends to a 
limit if xf \ 0. 

□ . 


Lemma A.4. As xf f tt, the function F(xf) tends to a limit. 

□ . The function F(xf) = F(xf), xf G [0,7r] where xf = n — xf fulfils identical equation like F(xf) that is why all 
qualitative properties of the function F(xf) pass to the function F(xf). In particular F(xf) tends to a limit as 
xf \ 0 what is equivalent to that F(xf) tends to a limit as xf | ir. 
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□ . 


Theorem A.5. As if \ 0 then every solution tends to a limit which may equal either kn or (k + \)n, where 
k £ Z. As if | 7 r then every solution tends to a limit which may equal In or (l + \)n, where l £ Z. 


Proof. Suppose that the limit lim^,_»o+ sin F(if), which exists by virtue of lem ma (A.3), is different from 
zero what means that sinF 0 7 ^ 0 where F a = lim^,^ 0 + F(if). From definition (Al) and from the lemma (A.2) 
we get then 


k F sin F sin if ~ sin F a + 4a, 


and further that 


8 a ■ 


k 2 sin 4 F n 


Ci 

if 2 ' 


2 k 2 sin 2 F 0 sin 2 if 

We must put 8 a + n 2 sin 4 F a = 0 otherwise F(if) would not be bounded because then F(if) = 0(\n.if) and 
it would contradict the lemma (A.3). But then F'(if) =o(if~ 1 ) and taking it into consideration in equation 
(A.2) multiplied by sin 2 if and taking the limit ^\0we get 


lim sin 2 ifF"(if) = sinF 0 cos F a . 

V>—>0+ 

Now if cos F 0 7 ^ 0 (sinF 0 7 ^ 0 from assumption) then F'(if) = 0(if _1 ) what cannot be because it contradicts 
just now derived formula F(if) =o(if~ 1 ). Thus if sinF 0 7 ^ 0 then cosF 0 = 0. So the conclusion is that the 
only possible limiting values of F(if) at zero may be looked for among the numbers from the countable set 
{fc-| : k £ Z}. As in previous proof we state that the behaviour of solutions in the left-sided neighbourhood 
of if = 7 r is similar i.e. the only possible limiting values of F(if) at n may be looked for among the numbers 
from the countable set {/■| : l £ Z}. This completes the proof. 

We are interested in behaviour of solutions of equation (A.2) in the vicinity of the singular points if = 0 
and if = n for which F(if) —■> kn with an integer k. Because of reflection symmetries of the equation we can 
limit ourself by considering only the solutions for which F(0) = 0 and F'( 0) > 0. For such solutions F(if) > 0 
and F 1 (if) > 0 for sufficiently small if. Irrespective of any behaviour of F as if —> 0 we can for contradiction 
take a sequence of points {^njneN for which if n —> 0 and F(if n ) > 0, F'(if n ) > 0. If for large enough n there 
was F(if) > 0 and F'(if) = 0, where if > if n , then from (A.2) we see that F"(if) would be positive and so F 
would have local minimum at this if. Thus for sufficiently small if we have always F(if) > 0 and F'(if) > 0. 
We shall use this fact fur ther . 

The formal expansion ©_ suggests that ifF'/F —> 1 as if —■> 0. This is the hint for the proof of asymptotics 
of solutions of the equation 


Lemma A. 6 . As if \ 0 then f(if) < 1 holds for the function f(if) := F ^ • 
□ . For contradiction suppose that F' sin 2if > 2 F then 


F' sin 2 if 
2 F 


F" sin 2 if F' , _ , ,. sin 2 if 

-1- - (2 F cos 2 if — F sin 2 if) < - —F, 

2 F 2 F 2 y ’2F 


since from assumption 2 F cos 2 if — F' sm2if < 2F(cos2if — 1) < 0. Moreover 

sin 2 F — (F') 2 sin 2 if < F 2 — \(F') 2 sin 2 2 if, 


so from (A.2) we get F" < 0. But then f(if) increases as if \ 0 and f(if) > c > 1 for if small enough. In 
particular F' > 2cF/ sin 2 if and integration in the interval (e n ,if n ) (where e n = if n /n) gives 


F(if n ) 2 F(S n ) 
> c- 


1 - 


F(e n ) 1 2 F(6 n ) 

- c - 

n sin 2 5 n 


1 - 


0 ^ d n <C ifm 


if n sin25„ \ nj e n n smzo„ \ 

and hence, taking the limit n —> 00 , we get g ^ eg for c > 1 where th e no nnegative number g is the limit 

g := lim^,_>o F('4 > )/' l l } - But this can be true only if g = 0 and then from (A.2) we get 

, sin 2 F f sin 2 F 


ifF" = 0(1) 


sin 2 F 


-0(l)F' + \ 


1 . 2 k 


if 2 if V sin 2 if 

fore some positive m. Integration of F"/F' < —m/if gives 


- (F') 2 < -mF' 


, F'(if) if 

fn-- < —mfn—, 

F'(e) e ’ 


0 < £ < if. 


But this implies that as £ \ 0 (at fixed if) then F'(s) diverges to +00 what in turn contradicts that 
F(if) =o(if). Thus the assumption F' sm.2if > 2F is false and always, for sufficiently small if, f(if) < 1 must 
hold. This proves the lemma. 
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□ . 


Lemma A.7. As ip \ 0 then f(ip) > 1 holds for the function f(ip) := 2 ^ ir [ . 

□ . For contradiction suppose that 2F' ta,nip < sin2F then 
' 2 F' tan ip 


sin 2 F 


2 tan ip 


2 F' 


(sin 2 F — F' cos 2 F sin 2 ip) > -^-t-Lf" 


sin 2F sin 2 2 F cos 2 ip K T ' sin 2 F 

since from assumption sin 2 F — F' cos 2F sin 2 ip > 2F' tan^>(l — cos 2 F cos 2 ip) > 0. Moreover 

sin 2 F — F' sin 2 ip > F'(2 tan ip — sin 2 ip) > 0, 
and sin 2 F — (F 1 ) 2 sin 2 ip > (F') 2 (4tan 2 ip — sin 2 ip) > 0, 


so from (A.2) we get F" > 0. But then f(ip) decreases as ip \ 0 and in fact f{ip) < c < 1 for ip small enough. 
In particular F' < c sin 2F/2 tamp and integration in the interval (E n ,ip n ) (where e n = ip„/n) gives 


F(ip n ) sin 2 F(5 n )ip n - e n 
< c- 


F(e n ) En ^ sin 2 F(S n ) F(s n ) 


< c 


2 tan Sn 


UEr, 


0 <C En ^ dn ^ Ipni 


ip n 2 tan 8 n ip n E n ip„ 

hence, taking the limit n —>■ oo, g ^ eg for 0 < c < 1 where the nonnegative number g is the limit g := 
lim^^o F(ip)/ip. But this can be true only if g = 0. The same arguments as in the previous lemma show that 
F'(ip) would be infinite in the limit ip \ 0 what would contradict that F{ip) =o(ip). Thus the assumption 
2 F' tan ip < sin 2 F is false and always, for sufficiently small ip, f(ip) > 1 must hold, what proves the lemma. 

□ . 

The inequalities we have already proved can be recast in to the form 

^ sin2jr w < ^F'M) < 

2 tan ip sin 2 ip 

and they are valid for sufficiently small ip. The Taylor series expansion about ip = 0 and F = 0 gives 

- -F 3 - -Ftp 2 +o(ip 3 ) < ipF' F < -iP 2 F +o{ip 3 ) 

3 3 3 

what implies that 

F{ip) = dip + 0(ip 3 ), where a = F'( 0), 

and this proves fj- The same of course, must hold at ip = n (the only difference is that F( n) = kn with 
an integer k, which is irrelevant to the proof). 
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